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A search for charmonium and other new states is performed in a study of exclusive initial-state- 
radiation production of DD events from electron-positron annihilations at a center-of-mass energy of 
10.58 GeV. The data sample corresponds to an integrated luminosity of 384 fb _1 and was recorded 
by the BABAR experiment at the PEP-II storage ring. The DD mass spectrum shows clear evidence 
of the i/>(3770) plus other structures near 3.9, 4.1, and 4.4_GeV/c 2 . No evidence for F(4260) -> DD 
is observed, leading to an upper limit of S(Y(4260) -> DD) / B(Y (4260) -> J/ipn+iT-) < 1.0 at 90% 
confidence level. 
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PACS numbers: 13.66.Bc, 13.87.Fh, 14.40.Gx 



The surprising discovery of new states decaying to 
has renewed interest in the held of 
charmonium spectroscopy, as the new states are not 
easy to accommodate in the quark model. In partic- 
ular, the BABAR experiment has discovered a new broad 
state, F(4260), decaying to J /il>n + ii~ in the initial-state- 
radiation (ISR) reaction e+e" -> 7/s_rF(4260). The 
quantum numbers J = 1 are inferred from the 
single virtual-photon production mechanism. Structure, 
possibly related to the y(4260), has been observed in 
the reaction e+e~ — ► "fisRip (2S)tt + tt~ 0]. A charmo- 
nium state at this mass would be expected to decay 
predominantly to DD, DD* or D*D* [1]. It is pecu- 
liar that the decay rate to the hidden charm final state 
J/ipir + ir~ is much larger for the Y~(4260) than for ex- 
cited charmonium states d, and that at the F(4260) 
mass the cross section for e + e~ — ► hadrons exhibits 
a local minimum Q. Many theoretical interpretations 
for the y(4260) have been proposed, including uncon- 
ventional scenarios: quark-antiquark gluon hybrids 0, 
tetraquarks @] and hadronic molecules [§[ . For a discus- 
sion and a list of references see, for example, Ref. (loj . 

This work explores ISR production of the DD final 
state for evidence of charmonium states and unconven- 
tional structures. A study by the BELLE collaboration of 
the DD*, and D*D* final states can be found in Ref. [11 1. 

This analysis is based on a 384 fb" 1 data sample 
recorded at the T(45) resonance and 40 MeV/c 2 be- 
low the resonance by the BABAR detector at the PEP- 
II asymmetric-energy e + e~ storage rings. The BABAR 



detector is described in detail elsewhere [12[. Charged 
particles are detected and their momenta measured by 
a combination of a cylindrical drift chamber (DCH) and 
a silicon vertex tracker (SVT), both operating within a 
1.5-T magnetic field of a superconducting solenoid. A 
ring-imaging Cherenkov detector (DIRC) combined with 
energy-loss measurements in the SVT and DCH are used 
to identify charged kaon and pion candidates. Photon 
energies are measured with a CsI(Tl) electromagnetic 
calorimeter (EMC). 

DD candidates are reconstructed in seven combina- 
tions of D decay modes, listed in Table [J [l3j]. In each 
channel we allow any number of photons in the event. 
Events are selected if the number of well-measured tracks 
is exactly equal to the total number of charged daugh- 
ter particles for the D and the D final states. Neutral 
pion candidates are formed from pairs of photons each 
having an energy greater than 30 MeV. The K° candi- 
dates are reconstructed in the tt + tt~ decay mode. The 
tracks of each D candidate are geometrically constrained 
to come from a common vertex. Additionally, for the 
D° — > K~ir + Tr° channel, 7r° and D° mass constraints are 
included in the fit, and for the D~ — > K®tt~ channel 



a Kg mass constraint is imposed. D candidates with a 
X 2 fit probability greater than 0.1% are retained. Sub- 
sequently, each DD pair is refitted to a common vertex 
with the constraint that they originate from the e + e~ in- 
teraction region; only candidates with a x 2 fit probability 
greater than 0.1% are retained. Extra tt° candidates may 
originate from random combinations of photons. Aside 
from 7r°'s from D° decays, we require that there be no 
more than one other ir° candidate in the event (except 
for channel 4, where we require that there are none). 

For D channels without a 7r° candidate, the D mo- 
mentum is determined from the summed 3-momenta of 
the decay particles and the energy is computed using the 
nominal D mass value 0, For the D° — > K~Tr + n° 
channel, the 4-momentum from the mass constrained fit 
is used. This procedure gives similar DD mass resolutions 
for all the channels. 

The ISR photon, preferentially emitted at small an- 
gles with respect to the beam axis, escapes detection in 
approximately 90% of events. We therefore reconstruct 
the ISR photon as a missing particle. We define the 
squared recoil mass (M 2 ec ) to the DD system using the 
four-momenta of the beam particles (p e ± ) and the recon- 
structed D (pd) and D (pjj): 



(Pe 



Pe 



Pd-Pd) 



(1) 



This quantity should peak near zero for ISR events and 
for exclusive production of e + e~ — * DD or e + e~ — * DD*. 
In the latter case, the DD mass distribution peaks at 
masses well above 6 GeV/c 2 . Therefore we select ISR 
events by requiring a DD invariant mass below 6 GeV/c 2 
and | M 2 ec |< 1 GeV 2 /c 4 . 

Monte Carlo simulations of e + e~ — ► jisrDD and can- 
didates from the process e + e~ — > ^isrJ/^P, J ftp - > 
K + K~tt + it~ in data are used to validate the require- 
ment on the number of residual tt° and the shape of the 
M 2 ec distribution. 

To estimate the number of background events in the 
signal region, the two-dimensional space spanned by the 
invariant masses of the two D candidates in each event 
is divided into nine regions: a central signal region and 
eight sideband regions above and below the signal re- 
gions, as illustrated in Fig. [T] for DD candidates recon- 
structed for the case of the K~n + and K + n~ modes. 
The mass range for the signal region is within ±2. 5a of 
the D mass, and the sideband regions are 2.5a wide and 
are separated from the signal region by 3.5a, where a 
is the mass resolution determined from a fit of a single 
Gaussian to the D candidate mass spectrum. 

The distribution of M 2 ec , summed over all DD chan- 
nels, is shown in Fig. O The shaded histogram corre- 
sponds to the background in the signal region estimated 
from the DD mass sidebands. The small inset in Fig. [2] 



5 



TABLE I: List of the reconstructed final states and corresponding values of efficiency times branching fraction. 



Channel First D decay mode Second D decay mode ef {m D p) (xlO 3 ) 



1. 


D°D° 


D° - 


-> K--K+ 




D° 


— ► K+-R- 


0.14 


2. 


D°D° 


D° - 


-> K~n+ 




D° 


— ► K+TT-TY 


0.42 


3. 


D°D° 


D° - 


-> K--K+ 




D° 


— > i\ + 7T~7r + 7r~ 


0.18 


4. 


D°D° 


D° - 


-> K--K+ 
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FIG. 1: _K" + 7r~ mass vs. 7f _ 7r + mass distribution for final 
state 1. The cross-hatched areas correspond to the signal and 
sideband regions. 



FIG. 2: Squared recoil mass, summed over all DD channels for 
ISR event candidates. The shaded histogram corresponds to 
non-DU background estimated from the DD-mass sidebands. 
The small inset shows the distribution of the center-of-mass 
polar angle of the DD system in the ISR region. 



shows the distribution of the DD center-of-mass polar an- 
gle 6 for DD candidates with | M 2 ec |< 1 GeV 2 /c 4 . The 
sharp peak at cosQ = — 1 is typical of ISR production 
and agrees with Monte Carlo simulations. 

The purity of each reconstructed D channel is demon- 
strated in Fig. [3] where projections of the candidate D 
mass distribution for events with | M 2 ec |< 1 GeV 2 /c 4 
are shown. Background is low in all channels. 

The DD mass spectrum summed over all channels (860 
events) is shown in Fig. [4] where the curves are the results 
from the fit described later. The shaded histogram repre- 
sents the background determined using the DD sideband 
regions and corresponds to 17.5% and 7.1% of the signal 
candidates for D°D° and D + D~ , respectively. We ob- 
serve a clear ?/;(3770) signal and other structures at the 
positions of i/'(4040) and i/>(4415). We also observe a sig- 
nificant structure in the 3.9 GeV/c 2 region, which may 
not be due to a resonance; the coupled- channel model 



of Ref. [15J in fact describes qualitatively the observed 
DD mass spectrum and the structure around 3.9 GeV/c 2 
without any need for additional ip states. 

To understand the background, we compute the ex- 
pected contribution from ISR production of the D*°D° 
system. Using Monte Carlo simulations and the cross 
section estimate from Ref. 11| we find ~ 6 % as possible 
contamination. This is confirmed by the examination of 
Z?7 and Dir° mass distributions where we find little ev- 
idence for D*° signal. In contrast, strong evidence for 
D*° production is observed for M 2 ec > 1.5 GeV 2 /c 4 . 
We investigate the possibility of background contribu- 
tions from DDX final states (where X ^ 7) by exploring 
events in the M 2 ec sideband region 1.5 < M 2 ec < 2.5 
GeV 2 /c 4 . The DD mass spectrum for these events shows 
no structure. We conclude that the residual background 
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FIG. 3: D-candidate mass projections for events with 
| M 2 ec |< 1 GeV 2 /c 4 and a DD invariant mass below 6 
GeV/c 2 . (a) K~ 7r + mass spectrum summed over channels 
1, 2, and 3. (b) K~tt + it mass spectrum summed over chan- 
nels 2 and 4. (c) _ft" _ 7r + 7r + 7r _ mass spectrum summed over 
channels 3 and 4. (d) K~7v + tv + mass spectrum for channel 5. 



to our signal is consistent with originating mostly from 
combinatorial non-DD events. 

In order to measure efficiency and DD mass resolution, 
ISR events are simulated at eight different values of the 
DD invariant mass between 3.75 and 7.25 GeV/c 2 . These 
events are generated using the GEANT4 detector simu- 
lation package [16} and are processed through the same 
reconstruction and analysis chain as are real events. The 
mass-dependent efficiency for each channel is fitted us- 
ing a second-order polynomial. The mass resolution is 
determined from the difference between generated and 
reconstructed DD mass. The DD mass resolution is sim- 
ilar for all channels and increases with DD mass from 
1.5 to 5 MeV/c 2 . We observe good agreement between 
Monte Carlo and data M^ ec distributions. 

We define Ni(m D jj) as the number of DD candidates 
for channel i. The channel branching fraction is Bt, and 
Ci(m D jj) is the efficiency as parametrized by the fitted 
polynomial. We define as ef(m D ^) the product efficiency 
times branching fraction for each channel, 
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FIG. 4: (a) The ISR DD mass spectrum. The shaded mass 
spectrum is from DD mass sidebands. The curve results from 
the fit described in the text, (b) An expanded view of the 
region with m D ^ < 4.2 GeV/c 2 . 



and then compute e (rn D jj) as 



( m z») = 



Ni(m L 



(3) 



The values of tf(ni D jj) are proportional to the expected 
yield for each channel. Their values, integrated over the 
DD mass spectrum, are reported in Table HI The result- 
ing yields, corrected for efficiency and branching frac- 
tions, are found to be consistent within the errors. 
The DD cross section is computed using 

dN/dm DD - 



^DD 



(in 



DD) 



(m D jj)dL I 'dm 



(4) 



DD 



7 



The differential luminosity is computed as 17 1 



dL 



dm 



= L 2m DD 



DD 



7TX 



-(ln(s/m 2 e )-l)(2-2x + x 2 ), (5) 



where a is the fine-structure constant, 



1 



DD 



/*, 



s is the square of the e + e~ center-of-mass energy, m e 
is the electron mass, and L is the integrated luminosity 
of 384 fb _1 . The background-subtracted cross sections 
for D°D° and D + D~ , averaged over 20 MeV/c 2 bins, are 
shown in Fig. [5J 
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FIG. 5: (a) D°D° and (b) D + D cross sections with statistical 
uncertainties only. 



Systematic errors on the cross sections (10.9% for 
D°D° and 8.1% for D + D~) include uncertainties in the 
particle identification efficiencies and tracking efficiency, 
possible inaccuracies in the simulation of extraneous 7r° 
candidates, and uncertainties in the background esti- 
mates (w 6 %) and on the luminosity function (w 1 %). 



Integrating the cross sections in the 0(3770) region 
(3.74-3.80 GeV/c 2 ), we compute the ratio of branching 
fractions, 



g(0(377Q) -> D°D°) 
£(V>(3770) -► L>+D-; 



= 1.78 ±0.33 ±0.24, (6) 



to be compared with the value of 1.28 ± 0.14 reported by 
the PDG Q. 

We perform an unbinned maximum likelihood fit to the 
DD mass spectrum summed over all channels. The pa- 
rameters of the V(4040), -0(4460), and -0(4415) are fixed 
to the values reported in Ref. [18| while the F(4260) 
parameters are taken from our measurement from the 
J/-07r + 7r~ channel 0. The parameters of the ip(3770) 
are left free in the fit. In addition, we search for evidence 
of the y(4260) in this spectrum. Resolutions effects have 
been ignored since the widths of the resonances are much 
larger than the experimental resolution. 

We express the total DD production as 



/ P + ciWie^ 1 + c 2 \/Ge 



i4>2 



■ + c n W n e i ^ +(l-f)B, 
(7) 

where Cj and <pi are free parameters, Wi are spin-1 rela- 
tivistic Breit-Wigner distributions, P represents the non- 
resonant contribution, B describes the non-DD back- 
ground and / (0.829 ± 0.015) is the signal fraction. The 
efficiency e B (m D j=j) is almost linear and increases from 
w 2 x 10~ 3 to w 4 x 10 -3 in the fitted mass region. It 
has been parametrized by a 2 nd order polynomial and 
it has been multiplied by P and Wi- The data require 
that wc include the 3.9 GeV/c 2 structure, as suggested 



Ref. 



15j, which we parameterize empirically as the 
square root of a Gaussian times a phase factor (VGe*). 
The parameters of the Gaussian are left free, and the 
phase allows interference with the tp states. 

We find that, in order to have a satisfactory description 
of the data, interference must be allowed between the 
resonances and the non-resonant contribution P. The 
latter contribution is parametrized either as a linear (a + 
bm) or a threshold function (to — m t f l ) a e~ bm ~ cm , where 
to = to D |j, m t h is the threshold DD mass, and a, b and 
c are free parameters. This threshold function has also 
been used to describe the non- DD background B. 

The two different parametrizations give similar results, 
which are considered in the evaluation of the systematic 
errors. These include also uncertainties in the D mass 
and on the overall DD mass scale. The size of the non- 
resonant production is determined by the fit. 

The fit with a linear non-resonant contribution is 
shown in Fig.[3](a). Figure HJb) shows an expanded view 
of the threshold region. 

The fit returns the following parameters for the 
G(3900) structure and for the ^(3770): 

to(G(3900)) - (3943 ± 17 stai ± I2 syst ) MeV/c 2 , (8) 



8 



a(G(3900)) = (52 ± 8 stai ± 7 syst ) MeV/c 2 , (9) 



m(V>(3770)) = (3778.8 ± 1.9 stat ± 0.9 syst ) MeV/c 2 , (10) 



r(^(3770)) = (23.5 ± 3.7 stat ± 0.9 syst ) MeV. 



(11) 



The systematic error on the 0j(377O) mass includes un- 
certainties in the D mass, background parametrization, 
and detector related issues such as magnetic field, EMC 
corrections and energy loss. We measure a significantly 
higher -0(3770) mass with respect to previous measure- 
ments (3772.4±1.1) MeV/c 2 [<S\. The change in likelihood 
due to the inclusion of a F(4260) amplitude in the fit is 
given by A(21n(i)) = 0.1 with two additional fit param- 
eters. 

The systematic errors due to the masses and the widths 
of the 0(4040), ^(4160), ^(4415) and y(4260) resonances 
in the fit are evaluated by varying them by their statis- 
tical uncertainties. The signal fraction has been varied 
within its statistical error and the meson radius used in 
the Blatt-Weisskopf damping factor [l|| present in the 
relativistic Breit-Wigner has been varied between and 
5 GeV -1 . The deviations from the central values are 
added in quadrature. The uncertainty on e B (m D ^) is 
evaluated by using a weighted mean of branching frac- 
tion and efficiency uncertainties for the different chan- 
nels. The fitted y(4260) yield before efficiency correction 
is 0.2 ± 6.1 Bt at ± 2.8 syst events^ 

This y(4260) yield in the DD channel is used to com- 
pute the cross section times branching fraction, which 
can then be compared to our measurement from the 
J/ipTr + n~ channel Q. We obtain 



6(F(4260) -> DD) 
£(F(4260) J/07T+7T-) 



< 1.0, 



(12) 



r(r(4260) 



e + e~ 



) ■ £(F(4260) -» DD) < 5.7 eV, (13) 



at 90% confidence level. 

In conclusion, we have studied the exclusive ISR pro- 
duction of the DD system. The mass spectrum is dom- 
inated by J PC = 1~~ states; in particular the ^(3770) 
is clearly seen. In order to fit the mass spectrum, sig- 
nals from -0(4040), -0(4160), and -0(4415) have been in- 
cluded. The fit requires the presence of a broad structure 
near 3900 MeV/c 2 . The presence of an enhancement in 
this region is predicted by a coupled channel model from 
Eichten et al. [15| . although the possibility of the pres- 
ence of a new -0 state cannot be excluded. 

If the F(4260) is a 1~~ charmonium state, it should 
decay predominantly to DD 4]; however no evidence is 
found for F(4260) decays to DD. Other explanations 
have been proposed, such as a hybrid, baryonium, or 
tetraquark state. 
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